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Abstract 

We compute the two- loop QCD helicity amplitudes for the process e + e~ — > qqg. The amplitudes 
are extracted in a scheme-independent manner from the coefficients appearing in the general tensorial 
structure for this process. The tensor coefficients are derived from the Feynman graph amplitudes by 
means of projectors, within the conventional dimensional regularization scheme. The actual calculation 
of the loop integrals is then performed by reducing all of them to a small set of known master integrals. 
The infrared pole structure of the renormalized helicity amplitudes agrees with the prediction made by 
Catani using an infrared factorization formula. We use this formula to structure our results for the finite 
part into terms arising from the expansion of the pole coefficients and a genuine finite remainder, which is 
independent of the scheme used to define the helicity amplitudes. The analytic result for the finite parts 
of the amplitudes is expressed in terms of one- and two-dimensional harmonic polylogarithms. 
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1 Introduction 



The three-jet production rate in electron-positron annihilation and related event shape observables were 
measured to a very high precision at LEP, where they were used in particular for the determination of the 
strong coupling constant a s . At present, the error on the extraction of a s from these data is dominated by 
the uncertainty inherent in the theoretical next-to-leading order (NLO) calculation ||-(7) of the jet observables 
(see H for a review). The planned TESLA 0] linear e + e~ collider will allow precision QCD studies at energies 
even higher than at LEP. Given the projected luminosity of TESLA, one again expects the experimental errors 
to be well below the uncertainty of the NLO calculation. 

The calculation of next-to- next-to-leading order (NNLO), i.e. G(a%), corrections to the three-jet rate in 
e + e~ annihilation has been considered as a highly important project for a long time ]T3). In terms of matrix 
elements, it requires the computation of three contributions: the tree level 7* — ► 5 partons amplitudes |ll| |l3[ |, 
the one- loop corrections to the 7*^4 partons amplitudes |l4|-|l7t , and the two-loop (as well as the one- loop 
times one-loop) corrections to the 7* — > 3 partons matrix elements. In a previous publication JL8||, we have 
derived both the interference of the tree and two- loop matrix elements and the self- interference of the one-loop 
amplitudes averaged over all external helicities. In the present work, we extend this calculation to compute 
the two-loop helicity amplitudes for the process e + e~ — > qqg. 

The most precisely measured observables related to e + e~ — > 3 jets are the jet production rate itself and 
a number of event-shape variables. The calculation of these phenomcnologically most relevant applications, 
which also dominate the extraction of a s , at NNLO accuracy requires only the helicity averaged squared 



matrix element at the two- loop level derived in 18 1. Nevertheless, the helicity amplitudes presented here are 
interesting for a number of reasons: 

• Oriented event-shape observables, which measure the spatial orientation of the final-state jets relative 
to the direction of the incoming beams require, even for unpolarized beams | fL9| , the calculation of the 
polarization tensor of the virtual photon mediating the interaction. This polarization tensor can be 
recovered from the helicity amplitudes. 



Likewise, to determine the direction of the decay leptons in the crossed process, V + 1 jet production 
at unpolarized hadron colliders, it is necessary to compute the polarization tensor of the vector boson. 

• Polarization of the beams is an important option for the future linear e + e~ collider TESLA 0, thus 
providing a direct measurement of event-shape observables in polarized e + e~ annihilation. 

• NNLO predictions for (V + l)-jet production at the RHIC polarized proton-proton collider and for 
(2 + 1) jet production at a currently discussed polarized upgrade of the HERA collider do require the 
calculation of the two- loop helicity amplitudes. These observables would then form part of a full NNLO 
determination of the polarized parton distribution functions in the proton. 

• The study of formal aspects of two- loop matrix elements, such as their collinear limits or their high 
energy behaviour can be carried out more elegantly on the basis of the underlying helicity amplitudes. 

Two-loop helicity amplitudes have up to now only been derived for 2^2 bosonic scattering processes 
with all external legs on-shell: for gg — > 77 ^0|, 77 — > 77 (2l],|2^] and gg — > gg p3|,p4|. The latter calculation 
also confirmed earlier results for the squared two-loop gg — > gg matrix element p5[ . In the above calculations, 
which were all carried out within dimensional regularization p6|-p8| with d — 4— 2e space-time dimensions, two 
different methods were used to access the helicity structure of the matrix element: explicit contraction with 
the external polarization vectors [^0|, ^l| , |23|,|2^1 or projection onto the individual components of the Lorentz- 
invariant decomposition of the amplitude fl22| . Once these are applied to expose the helicity structure, one is 
left with the task of computing a large number of two-loop integrals. Using integration-by-parts [p8|-[3C|l and 
Lorentz-invariance |u] identities, these can be reduced |3^] to a small number of so-called master integrals, 
which were derived for massless on-shell two- loop four-point functions in If an explicit contraction 

with the external polarization vectors is performed, one also has to compute two-loop integrals over the (d— 4) 
dimensional subspace of loop momenta, which reduce however to simple vacuum diagrams . For 2 — > 2 
scattering processes with external fermions and all external legs on-shell (e + e~ — > e + e~ , qq — > q'q' , qq — > qq, 
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qq — > gg_,qq — ^ g"f and qq — > 77), only the squared, helicity-averaged two- loop matrix elements were computed 
so far |39H42[|- 

The method employed here to extract the two-loop helicity amplitudes for e + e~ — > qqg is similar to the 
approach of p2| by applying projections on all components of the Lorentz-invariant decomposition of the 
amplitude. Using this approach, the corresponding one- loop helicity amplitudes were derived in ||. The 
master integrals relevant in the present context are massless four-point functions with three legs on-shell and 
one leg off-shell. The complete set of these two-loop integrals was computed in jl3| , while earlier partial results 
had been presented in The master integrals in are expressed in terms of two-dimensional harmonic 

polylogarithms (2dHPLs). The 2dHPLs are an extension of the harmonic polylogarithms (HPLs) of jf8|. All 
HPLs and 2dHPLs that appear in the divergent parts of the planar master integrals have weight < 3 and can 
be related to the more commonly known Nielsen generalized polylogarithms p9|,|50|] of suitable arguments. 
The functions of weight 4 appearing in the finite parts of the master integrals can all be represented, by 
their very definition, as one-dimensional integrals over 2dHPLs of weight 3, hence of Nielsen's generalized 
polylogarithms of suitable arguments according to the above remark. A table with all relations is included 
in the appendix of [Q. Numerical routines providing an evaluation of the HPLs and 2dHPLs [^2| are 
available. 

After carrying out ultraviolet renormalization of the amplitudes in the MS scheme, one is left with poles 
which are purely of infrared origin. The infrared pole structure of the amplitudes can be predicted using 
Catani's infrared factorization formula psj . We use this formalism to present the infrared poles and the finite 
parts of the helicity amplitudes in a compact form. 

This paper is structured as follows: in Section ^, we outline the calculational method used to derive the 
helicity amplitudes and discuss the techniques used to extract the ultraviolet and infrared pole structure. 
We also elaborate on the relation to previous work. The two-loop helicity amplitudes are computed (in the 
Weyl-van der Waerden formalism, which is briefly described in the Appendix) in Section [|. Finally, Section ^] 
contains a discussion of the results and conclusions. 

2 Method 
2.1 Notation 

We consider the production of a quark-antiquark-gluon system in electron-positron annihilation, 

e + (f5) +e~(p 6 ) -> 7*(p 4 ) — ► q(pi) + qipi) + .9(^3) • (2.1) 
It is convenient to define the invariants 

S12 = (pi +P2) 2 , S13 = (pi +P3) 2 , s 23 = (p 2 + P3) 2 , (2.2) 

which fulfil 

Pi = (Pi +P2+ P3) 2 = Sl2 + S%3 + S 2 3 = S123 , (2.3) 

as well as the dimensionless invariants 

x = S12/S123 , y = S13/S123 , z = s 23 /si23 , (2.4) 

which satisfy x + y + z = 1 . 

The renormalized amplitude \M) can be written as 

\M) = V^S^q-g-q) , (2.5) 

1 Notc that an alternative approach avoiding the need to use the integration-by-parts and Lorentz-invariance identities to 
reduce the integrals appearing in the Feynman diagrams to a basis set has recently been proposed |46l (47| . This method relies on 
obtaining analytic expressions for the basic topologies with arbitrary powers of the propagators andarbitrary dimensions, which 
can often be found in terms of nested sums involving T-functions. The T-functions can be directly expanded in e and the nested 
sums related to multiple polylogarithms. 
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where V^ 1 represents the lepton current and denotes the hadron current. In a previous paper |18| , we have 
considered the unpolarized decay process 

1*(P4) — » q(pi) + qfa) + 9iP3) • (2.6) 



for which the amplitude is obtained from Eq. (2.5) by replacing the lepton current by the polarization vector 
of the virtual photon e^. The hadron current may be perturbatively decomposed as 

S,(q; g; q) = VW 9 v^I T% Uf (g; g; q) + (g) (q; g; q) + (g)' (q; g; q) + 0(a 3 s )\ , (2.7) 

where e q denotes the quark charge, a is the adjoint colour index for the gluon and i and j are the colour 
indices for quark and antiquark. a s is the QCD coupling constant at the renormalization scale /i, and the 
are the i-loop contributions to the renormalized amplitude. Renormalization of ultraviolet divergences is 
performed in the MS scheme. 

2.2 The general tensor 

The most general tensor structure for the hadron current 5 M (g; g\ q) is 

S^(q; g; q) = u{pi)i> 3 u(p 2 ) (Ane 3 .pi p lf2 + A 12 t 3 .p x p 2fl + Ai 3 e 3 .pi P3^) 

+ u{p 1 )i> 3 u(p 2 ) {A 21 e 3 .p 2 pijj, + A 22 e 3 .p 2 p 2fJ , + A 23 e 3 .p 2 p 3/x ) 

+ u{p\)l^u{p 2 ) (Bie 3 .pi + B 2 e 3 .p 2 ) 

+ u{pi)i 3 u{p 2 ) (Cipifi + C 2 p 2 ^ + C 3 p 3M ) 

+ DiU (pi)^37^(P2) 

+ Dsuip^fofaufa) , (2.8) 

where the constraint e 3 ■ p 3 = has been applied. All coefficients are functions of si 3 , s 23 and S123. The 
above tensor structure is a priori <i-dimensional, since the dimensionality of the external states has not yet 
been specified. The hadron current is conserved and satisfies 

S»{q;g;q)p% = 0; (2.9) 

it must also obey the QCD Ward identity when the gluon polarization vector €3 is replaced with the gluon 
momentum, 

S^(q;g;q)(e 3 ^p 3 ) = 0. (2.10) 

These constraints yield relations amongst the 13 distinct tensor structures and applying these identities gives 
the gauge-invariant form of the tensor, 

5 M (g;g; q) = ^4n(si 3 , s 23 , si 23 )Tn^ + Ai 2 (si 3 , s 23 , si 23 )Ti 2fl + Ai 3 (s 13 , s 23 , s 123 )T 13lJ , 

+ A 21 (si 3 , S 23 ,Si 23 )T 2 i^ + A 22 (si 3 , S23, 8123)^22^ + ^23(S13, S23, Si23)T 23M 

+ B(s 13 ,s 23 ,s 123 )T^ (2.11) 
where Aij and B are gauge-independent functions and the tensor structures T/j^ and T M are given by 

Tijf, = u(Pi)h u iP2)£3-PiPj t , - ^Yu(pi)<f. 3 u{p 2 )pj^ + ^-u(p 1 )fop' 3 j ll u(p 2 ) 7 (2.12) 

T 2 .jn = u(px)fcu(p 2 )es-P2Pj» - -Y-u(pi)fou(p2)pj» + -^w(pi)7^ 3 / 3 w(p2), (2.13) 

1 

2' 
1 

2'' 



Tft = s 23 ^u(pi)7 A1 u(p 2 )£3-Pi + 7; u (Pi)fafcl^u(p 2 ) 

- s 13 ( u(pi)7 At u(p 2 )e3-P2 + ^m(pi)7m^3^m(p2) ) ■ (2.14) 
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Each of the tensor structures satisfies both current conservation and the QCD Ward identity. The coefficients 
are further related by symmetry under the interchange of the quark and antiquark, 



^2l(si 3 ,S23,Si2 3 ) = -^12(523,513,5123), 

-A22(S13, S 23 , S123) = -^4ll(S23, 513,5123), 

^23(si3,S23,Si2 3 ) = -^13(S23, 513, 5123), 

-B(Si3,S 2 3,Si2 3 ) = B(s 2 3,S 13 ,S 123 ). (2.15) 

2.3 Projectors for the tensor coefficients 

The coefficients Ajj and B may be easily extracted from a Feynman diagram calculation, using projectors 
such that 

E^ffl e^(q;g;q)=X(s 13 , s 23 , 5i 23 ). (2.16) 

spins 

The explicit forms for the seven projectors in d space-time dimensions are, 

V(A \ (523Sl23rf+Si 3 Sl2(d- 2)) f „ (s 13 + s 23 ) (d - 2) t # 

" 2 S 3 3S 2 2 ( d _3) Sl2 3 in ' £4 2^ 3S 2 2(d _ 3)si23 J 12^4 

((S23 + 5l2)d + 2si 3 ) | (s 2 35l2 3 (rf- 2) + Si 3 Si 2 (d - 4)) ^! 



2 Sl2S ? 3 (d-3) Sl2 3 13 4 2 S 235? 2 5? 3 5l23(d-3) " 21 

(5l3 + 5 23 )(rf-4) rf (523 + Sl2)(rf-4) ^ 

2(rf-3)s 2 2 Si23Si 3 S23 22 4 2 S23 S 12 S 2 3 S 123 (d-3) 23 4 
1 



25 2 3 5 2 2 (d-3) J ■ £4 ' 

pM x (ai 3 + a 23 )(d-2) f . (d - 2)( S23 si 2 (rf - 4) + sj^jd - 2)) f , 

{Al2) ~ 2 S 2 3S f 2 (d-3) Sl23 11 4 2 S 2 3S 2 2S2 3(rf-3) S i23(d-4) Jl2 ^ 4 

(d-2)(s 13 + s 12 ) t * , ((rf-6)(rf-2)( 3l3 + S23)-43i 2 ) t , 

25 2 3Sl2 S23(d-3) Sl2 3 13 " 4 2(d-4) S 2 2Sl 35235 12 3(d-3) 21 ' 4 

(523512(^-4) + s 13 s 123 (d- 2)) f , (2s 23 + (ai3 + si 2 )(d-2)) t , 
25 2 2S i 3 5 2 3 (d-3) Sl23 J22 ' £ 4+ 2 Sl2Sl3S 2 3S i 23 (d-3) 123 



(d-2) 



-Tt- 



c 4- 



2(rf-4) S 2 2 si3S 23 (d-3)" 
((s 2 3 + si2)d + 2si 3 ) rt g „ _ (rf-2)(ai 3 + si 2 ) r t 
2s 12 s 3 3 (d-3)si23 U ' 64 2s 2 3 si2S23(d-3) Sl 23 12 



(ai 3 a 2 3(rf-2) + ai 2 ai 23 d) r t £ * ((fig + 5 23 )(d - 2) + 2s 13 ) yf 
2sf 3 si2S 23 si 2 3(c?- 3) 13 4 2s 2 3 si 2 s 23 (d- 3)s i23 21 

(513 + 512) 4) rpt * (513 + Sl 2 )(5 2 3 + 5l 2 )(rf - 4) t „ 



' 2si2Si3S 2 3 si23(d- 3) 22 4 ~ 2s 2 3 si 2 si 3 si23(d-3) 23 
j \ t4 . f * 

2 S2 35 2 3 5l2(d-3) 4 ' 

pM v (g23ai 23 (d-2) + 8i 3 si 2 (d-4)) t , (-4ai 2 + (ai 3 + s 23 )(d-6)(rf-2)) t , 

1 21J ~ 2 S 2 3S 2 2 s 23 (d-3) Sl23 J "- £ 4+ 2(d-4) S 2 2Sl3S23Sl23 (d-3) Jl2 ^ 4 



| (s 23 + 5i 2 )(d-2) + 2s 13 ) rt e . | (d-2)(s 2 35i23(rf-2) + 5i35i 2 (d-4)) f 
2s 2 3 si2S23(d-3)si2 3 13 ' 64 2sf 2 Si3s| 3 si23(d-3)((i-4) 
(5i 3 + 5 23 )(rf-2) rt (g 23 + ai2)(d-2) Tt 

2s 2 2 s| 3 si23(rf-3) 22 ' £4 2si3si 3 si2(d-3)si23 23 ' 64 

1 ( rf - 2 ) rt . e * 



2(rf-4)s 2 2 si 3 S23(d-3) 
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P(A 22 ) 



(si 3 +S 23 )(d-4) Tt ^ 
2S13S?2 S 23(^- 3)S123 11 ! 



«4 - 



(s23Sl2(rf- 4) + Sl 3 Sl23(rf- 2)), 



2sj 2 s 13 sl 3 (d- 3)si2 3 

0>13 + ai2)(rf-4) t » _ (£i3+£23Kd-2) t „ 
2 D M J 13' e 4 r> o 2 „ oN J 21' e 4 



12 c 4 



+ 



2si 3 s^ 3 si 2 (d- 3)si2 3 

(S23Sl2(d- 2) + Si 3 Si2 3 rf) t 



2s 12 S 23 S 123(^ - 3)" 



T(A 23 ) = 



2s\ 3 s\ 2 {d- 3)si2 3 
2 1 T"t . e * ; 

2S23 S 12 — 3) 
(s 23 + Sl2)(rf-4) ^ 



(si 3 d + si 2 d + 2s 23 ) r 
2si 2 s% 3 si 23 (d - 3) 



- 23 ' e 4 



e 4 + 



{2s 23 + (s 13 + s 12 )(d - 2)) ^ 



2s 23 si2sf 3 si 23 (d- 3)" 11 " 4 ' 2si 3 sl 3 s 12 (d - 3)s 123 

(S13 + 512)(523 + Sl2)(d - 4) yt t 

2sj 3 s 12 s^ 3 s 123 (d - 3) 13 4 2si 3 s^ 3 si 2 (d - 3)si 23 
((£13 + si 2 )d + 2s 23 ) T t * , (si 3 s 23 (d - 2) + s 12 si 23 d) 
2s 12 s 3 23 s 123 (d-3) 22 '' 4 



1 12 e 4 



g 23 +5l2)(rf-2) | # 
'21 ' f l 



Try * i V wio v~ — / i " i-zo ~y rp] * 



1 



yt . e * 

2s^ 3 si 3 si 2 (d- 3) 4 ' 



P(B) 



2(d - 3)s 2 2 s 2 3 
(d-2) 



1 ll e 4 



2si 3 s 12 s3 3Sl 2 3 (d-3) " 23 4 
(d-2) 



T t 

2(d-4)s 2 2 s 13 s 23 (d-3) 12 

1 



e 4 + 



2s 23 si 2 (d- 3)sf 3 



J 13 e 4 



2(d-4)sf 2Sl3S23 (d-3) 

1 yt . e * 

2(d- 4)s 2 2 si 3 s 23 



T 2 \ ■ e\ 



2 s 2 3 S 12(d 3) 



J 22 £ 4 



2si 3 si 2 (d - 3)s 2 3 



2 '23 c 4 



(2.17) 



2.4 The perturbative expansion of the tensor coefficients 

Each of the unrenormalized coefficients Ajj and B has a perturbative expansion of the form 

2 



A}j = VlWqVi^I T" 



,(0),un 



(£)43" + (s) 

B (0),un + / ^ B (l),un + / B (2),un + ^3) 



where the dependence on (si 3 , «2 3 , S123) is implicit. At tree level, 

^°J' un (s 13 ,S2 3 ,Si2 3 ) = 0. 

B (0) ' Un (si 3 ,S2 3 ,Sl2 3 ) = 



Sl 3 S2 3 



(2.18) 

(2.19) 
(2.20) 



The one-loop contributions can be written in terms of the one-loop box integral in d = 6 — 2e dimensions, 
Box 6 (sjj, Sik, Sijk), and the one-loop bubble, Bub(sy ), as follows: 



,(l),un/ \ 

(S13, «2 3 , S123) = 

(d-4) 



AM 



Bub(si 23 ) - [ d ^ [Bub(si 3 ) - Bub(si 23 )] 



2(si 3 + si 2 )si 3 v 2si 2 si 3 

((d - 2)s 23 si2 + (d- 4)s 23 si 3 + 4si 2 (si2 + s 13 )) 
2si 2 si 3 (si3 + S12) 2 

Box 6 (s 13 ,s 23 ,si 23 ) 



[Bub(s 23 ) - Bub(si 23 )] 



(d-4)(4si 2 + (d-2) S23 ) . 



4si 2 si 3 
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. (si2 + si 3 )(ds 2 3 + 4s i3 ) + 2s 23 si3 



2(S1 3 + Si 2 ) 2 S? 3 



[Bub(s 23 ) - Bub(si 23 )] 



+ (d "^3~ 6) Bax»(- ia , - 18) «i M ) + (^Kg|±i£H) B ox«( Sl2 , S23 , Sl23 ) j , (2.21) 

^12 ),Un ( s 13,S23,Si 23 ) = 

ivf- {( j~ W) , [Bub(s 13 ) - Bub( Sl23 )] - (( rf ~ 10) , SlS ~ 4S1 ; } [Bub( S23 )-Bub( Sl23 )] 

+ -10)-) B„x»( SI 3, 32s , sl23) ) 

4si 2 si 3 y 

+ 7; ; ; r Bub(s 23 ) - Bub(s i23 ) 

2si 2 si 3 (si 3 + s 12 ) 

((d - 2) 2 s 12 s 13 + 4(d - 4)s 12 s 23 + 2(d - 4)(d - 6)s 13 s 23 ) 6 

H : BOX (Si2,Si3,Si23) 

4s i2 si 3 s 23 

4si 2 si 3 J 

^13 ),Un ( s 13,S23,Sl23) = 

N ( 2(si f + "^ )si3 [Bub( S23) - B U b( Sl23 )] + (<f -^- 6) Box«( Sl3 , S23 , Sl23 )) 
1 '" ^ ! Bub( Sl23 ) - [Bub(s 13 ) - Bub( Sl23 )] 



N\ si 3 (s 23 + si 3 ) 2s 23 si 3 

(4si 2 s 2 3 - dsi 2 (s i3 + S23) 2 - 2(d - 2) si 3 s 23 (5i 3 + s 23 )) 

2(s 23 + S13) 2 S 2 3 S23 



+ ^12^13 - ^12^13 -r 02 3, ~ ^13^23^13 X *2 3 ;; [Bub(si2 ) _ Bub(si23)] 



(2(d- 3)313 + tfa 12 ) [R , , , p, . , (d-4)((d-2) Sl2 +4s 23 ) 6 

—2— ■ — Bub(s 23 ) -Bub(si 23 ) Box (si2,si3,si23) 

2s i 3 ( s i3 + S12) 4s 23 si 3 

(d-2)(d Sl2 + 2(d-J) £l 3) Box6(si2)S235 Si23) \ ( ^ 3) 



4s 2 



S (lKUn (si 3 , S23 , S l 23 ) 



/ d 2 -3d + 4 

" V-^ E " b(s "'» 

+ (4(d-3)y i (,, a + ^) + W-4)( d -7 )w ,3) p ub( , ia) _ Bub(si23)] 
2si 2 s 23 (s 23 + si 2 )si 3 (d - 4) 

+ (4(d-3).,f 2 + ( d -2)(d-7) w ,,) Box « (si3> , 23 , ai23) ^) 

»Sl2Sl 3 S 23 J 

1 /(7d- 16- d 2 ) „ , , , (16 -5d) r „ , . , „ , , 
+ N 4(d-4) Sl3 J + 4 ( d-4) Sl3 ; 23 [Bub(Sl2) - Bub(Sl23)] 
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(gig + (d- 6)s 23 ) 

2(s 23 + S 12 )s 2 3Si2 



[Bub(s 13 ) - Bub(si 23 )] 



(4(d - 3)s 23 si2 + (d - 4)(d - 6)si 3 s 23 + (d - 2)s 12 s 13 ) 6 



4si 2 si 3 s 23 



Box (S12, S13, S123) 



S13 <-> S 23 



(2.24) 



Explicit expansions of the one-loop integrals around e ~ in terms of HPLs and 2dHPLs are listed in 
Appendix A of g§. 

Similarly, the unrenormalized two-loop A 2 j' nn and B^ ,un coefficients were obtained analytically (making 
extensive use of the computer algebra programs MAPLE Q, FORM2 Q and FORM3 Q, where the latter 
two are particularly well suited for handling the large-size expressions arising at intermediate stages of the 
calculation) in terms of a basis set of two-loop master integrals. This basis set comprises 14 planar topologies 
and 5 non-planar topologies. Five of the topologies require more than one master integral, so that in total 
24 master integrals are needed. A more detailed discussion can be found in Ref. p8[ . However, we note that 
Laurent expansions for each of these master integrals have been derived in |43| ] by solving differential equations 
for the master integrals (equations that are differential with respect to the momentum scales involved in the 
diagram). The e-expansions of 4 2 j' Un and B( 2 )' un can therefore be obtained by directly substituting the 
e-expansions of the individual master integrals. 



2.5 Relation to previous work 

We have considered the case where the correlations with the lepton current are ignored in a previous paper Jlq ] . 
In this instance, the squared amplitude for the process 7* — > qqg, summed over spins, colours and quark 
flavours, was denoted by 

(M\M) = ]T |e 4 • S(q;g;q)\ 2 = T(x,y, z) . (2.25) 
The perturbative expansion of T(x,y, z) at renormalization scale ji 2 = q 2 = S123 reads: 



T{x,y,z) = 16TT 2 a^2e 2 a s (q 2 



^l^T^ {x ,y, z) + 0(al(q 2 )) 



where 



T^ 2 \x,y,z) = (M W \M {0) ) =4V(l-e) 



(1 _ e ) 1 y _|_ i) + 2 ^ z y z z ) ~ 2eyz 

z y) yz 



T^(x,y,z) = (M {0) \M W ) + (M (1) \M (0) ) , 

T^(x,y,z) = (M W \M W ) + (M (0) \M {2) ) + (M^\M W ) 



(2.26) 



(2.27) 

(2.28) 
(2.29) 



where V = N 2 — 1, with N the number of colours. T^ 4 \x,y,z) was first derived in HI through to O(e ) 
while an explicit expression for it to all orders in e was given in [^8[ . The contribution to T^(x, y, z) from 
the interference of two-loop and tree diagrams 



T^ 2x °]}(x,2,,z) = (M^\M^) + {MW\M<V) , 
as well as the one-loop self-interference 

T^ 1 ^\x,y 1 z) = (M^\M^) 

were first derived in 11181. 



(2.30) 



(2.31) 
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It is straightforward to obtain the interference of the tree and i-loop amplitudes in terms of the tensor 
coefficients, Ajj and B. We find 

{M {0) \M {1) ) = 

^2(1 - e) ((si 2 si23 + S12S13 + S13S23) - e(si3 + s 2 3)(si 2 + S13)) A^(si3, S23, S123) 
+ (2(si2 + S23) 2 - 2e (S123S23 + (S12 + S23) 2 ) + 2e 2 (si 3 + s 23 )(si2 + S23)) A$ (s 13 , s 23 , s 123 ) 

+ 2 (s 2 3 - e(si 3 + S 2 3))(si23 - e(si3 + S23)) A§(si3,S23j S123) 

+ 2 (s 13 + S23 + 2si 2 si23 - 2e(s 123 - S12S13 - S12S23 - S13S23) + e 2 (si 3 + S23) 2 ) -B W (si3, s 2 3, S123) 

+{pi«-> (2.32) 

The above relation holds for the unrenormalized as well as for the renormalized matrix element, involving 
the appropriate unrenormalized or renormalized tensor coefficients respectively. Similar, but more lengthy, 
expressions can easily be obtained for the interference of i- and j-loop amplitudes. We have checked that 
inserting the expressions for Aj l j and flW into Eq. (2.32) reproduces our earlier results [p| at the one- and 



two-loop level both at the master integral level and after making an expansion in e. 
2.6 Ultraviolet renormalization 

The renormalization of the matrix element is carried out by replacing the bare coupling ao with the renor- 
malized coupling a s = a s (^ 2 ), evaluated at the renormalization scale fj, 2 



where 



-T©+(f-§)(i;r+^> 



(2.33) 



S e = (47r) £ e" e7 with Euler constant 7 = 0.5772 . . . 

and /jIq is the mass parameter introduced in dimensional regularization p6|~p8| to maintain a dimensionless 
coupling in the bare QCD Lagrangian density; (3o and (3\ are the first two coefficients of the QCD /3-function: 

. 11C A -AT R N F 17C 2 a -10C a T r N f -6C f T r N f 

P0 — , Pi = 7, , ( z -^ 

6 6 

with the QCD colour factors 

iV 2 - 1 1 

C A = N, C F = ^ r , Tr =2- (2 ' 35) 

We denote the i-loop contribution to the unrenormalized coefficients by A^j' un and £?W- un , using the same 



normalization as for the decomposition of the renormalized amplitude (2.7); the dependence on (S13, S23, S123) 



is always understood implicitly. The renormalized coefficients are then obtained as 



(0) 


= 0, 




/./ 




(1) 

I.J 


= s; 


1 4(1), no 


(2) 
/,/ 


= s; 


2 4 (2),un 



2e 

and 

5(0) _ B (0),un 



A\ j = S- z AY>' uu -^Sr L AYj' uu , (2.36) 
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B (l) 



g-l£j(l),un _ PO g(0),un 



D 



(2) _ 



o-2 R (2).un 3/3 „_i R (l),un 3 ^o\ R (0). 

^ ^ _ "2T Ae ~ ^47" 8e^ 



(2.37) 



For the remainder of this paper we will set the renormalization scale /i 2 = q 2 . The full scale dependence 
of the tensor coefficients is given by 



Au = 
B = 



q s (a< 2 ) 



.4 



(i) 



2tt 



« s (m 2 ) 

2tt 



2 r 



2 r 



flW + f 3^ B (D + ^ B (0) ) In f ^ U In 2 f £ 



v47rae, 

2.7 Infrared behaviour of the tensor coefficients 



q 2 



3(3 2 



+ 0(ai)}. (2.38) 



After performing ultraviolet renormalization, the amplitudes still contain singularities, which are of infrared 
origin and will be analytically cancelled by those occurring in radiative processes of the same order. Catani [fj3[l 
has shown how to organize the infrared pole structure of the one- and two-loop contributions renormalized 
in the MS scheme in terms of the tree and renormalized one-loop amplitudes. The same procedure applies to 
the tensor coefficients. In particular, the infrared behaviour of the one-loop coefficients is given by 



A {1) 



A 



(1), finite 
IJ ' 



j(l)( e ) jB (0) +B (l),finite j 



(2.39) 



while the two-loop singularity structure is 



.4 



(2) 



lM(e)A?]+A 



j (2), finite 



B W = [-ij(i)( e )j(i)( e )-^J«( e ) 



_ £ ^r(l-2e) f(3 



+lW( f )B< 1 '+B (2)4 ' litc , 
where the constant K is 



K 



67 
18 



TV" 



r(i 



C A 



— + K) (2e) + (e) B 



(0) 



10. 
~9 



T R N F . 



(2.40) 
(2.41) 



The finite remainders A 



(i) , finite 
IJ 



and ijW> nnltc remain to be calculated. 



For this particular process, there is only one colour structure present at tree level which, in terms of the 
gluon colour a and the quark and antiquark colours i and j, is simply T°-. Adding higher loops does not 
introduce additional colour structures, and the amplitudes are therefore vectors in a one-dimensional space. 
Similarly, the infrared singularity operator I^\e) is a 1 x 1 matrix in the colour space and is given by 



* (1) (e) 



e7 



2r(l-e) 
where (since we have set /i 2 = S123) 



N 



3 



A) 
2Ne 



1 

N 



3 

27 



S123 



(2.42) 



(2.43) 



Note that on expanding Sjj, imaginary parts are generated , the sign of which is fixed by the small imaginary 
part +i0 of Sij. The origin of the various terms in Eq. ( |2.42j ) is straightforward. Each parton pair ij in 
the event forms a radiating antenna of scale s^. Terms proportional to S^- are cancelled by real radiation 
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emitted from leg i and absorbed by leg j. The soft singularities 0(l/e 2 ) are independent of the identity of 
the participating partons and are universal. However, the collinear singularities depend on the identities of 
the participating partons. 

Finally, the term of Eq. (2.40) that involves H^ 2 \e) produces only a single pole in e and is given by 



e>«7 



4eT(l 



(2.44) 



where the constant is renormalization-scheme-dependent. As with the single-pole parts of J^(e), the 
process-dependent can be constructed by counting the number of radiating partons present in the event. 
In our case, there is a quark-antiquark pair and a gluon present in the final state, so that 



H m = 2H {2) 



where, in the MS scheme: 



q 



5 


ll7T 2 \ 


"12 + 


144 ) 


409 


Uir 2 


" 864 ~ 


96 


25 




~ 216 





27^ +, -=- 



so that 



The factors H\ and Hg A> are directly related to those found in gluon-gluon scattering 



7T 

48 



4C3+ 432 



N + 
l)N F 



a 



72 
41 
108 



89 
108 

2 



NNp 



7T 

96 



N F 

V A 
2^ 3 32 



1 

iV2 



N 



(2.45) 



(2.46) 



(2.47) 



11tt z 
^72" 



N 



41 

54 



T 



i 



1 NNp + ) — — H 7V|. 

18 36 J V 54 24 / TV 27 F 



(2.48) 



r(2) 



quark-quark 



scattering |40| and quark-gluon scattering J4l| (which each involve four partons) as well as in the quark form 
factor |57]-|6C[| and gluon form factor |61jj . We also note that (on purely dimensional grounds) one might 
expect terms of the type S 2 j to be present in H^ 2 \ Of course such terms are 1 + 0(e) and therefore leave the 
pole part unchanged, only modifying the finite remainder. At present it is not known how to systematically 
include these effects. 



3 Helicity amplitudes 

We can extend the results of the previous section to include Z boson exchange, 

e + (P5) +e~(Pe) -> (Z*,l*) ipi) — * q(Pi) + q(P2) + g(ps) , 



(3.1) 



where the off-shell vector boson now distinguishes between left- and right-handed fermions by keeping track 
of the helicity of the final state quarks. 2 A convenient method to evaluate the helicity amplitudes is in terms 
of Weyl-van der Waerden spinors, which is described briefly in Appendix ^ and in detail in [p2||63[. 

It is also straightforward to include the spin-correlations with the initial state by contracting the hadronic 
current with the lepton current for fixed helicities of the initial state electron (and positron) . Using the 
spinor calculus of Appendix [A| we can express the lepton current in terms of the helicities of the incident 
e + and e~ (with momenta p§ and p$ respectively). Explicitly, 



K7(e++,e--) = e, 



L1„ 



V P&AP5B- 



V?(e++,e--) = e*j B 



Li 



P 6 AP5B- 



M 2 + iT z M z : 



Note that the full matrix element for any process should be summed over both photon and Z-boson exchange. 
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A R 

ec V p 5 apgb 



V?(e+-,e +) = ea* B p 5AP61 



s - M 2 7 + iV z Mz ' 



(3.2) 



The hadronic current is related to the fixed helicity currents, S Ab , by 

S^q+;g\;q-) = R^ ih V2aj B S AB (q+ ] g\; q-), 
(q-;gX;q+) = L\ h y/2 <yj B S AB (q- ;gX;q+). 



(3.3) 
(3.4) 

As in Eq. (2/7), the gauge boson coupling is extracted from S Ab . As mentioned earlier, the left- and right- 
handed currents couple with a different strength when the vector boson is a Z. 
The currents with the quark helicities flipped follow from parity conservation: 



Sab(i-\9^,Q+) = (S BA (q+;g{-xy.,q-))* . 

Charge conjugation implies the following relations between currents with different helicities: 

s AB(l X q''9^Q^q) = (-l)S AB (qXq; gX; qX q ). 

All helicity amplitudes are therefore related to the amplitudes with X q = + and = — . 
Explicitly, we find 



(3.5) 
(3.6) 



5 As(?+;fl'+;?-) = a(y,z) 



P1ADP2 P2B 

(P1P3) (P3P2) 
(pm)* 



(PlP3)(PlP2) : 



.^ z) ^ D P^b P^bP^A 

{PlP3){P3P2) KP1P3) KP3P2)* 

iPiAB + P2AB + PsAb) ■ ( 3 - 7 ) 



The other helicity amplitudes are obtained from S AB (q+; g+; q—) by the above parity and charge conjugation 
relations, while the coefficients a, /3 and 7 are written in terms of the tensor coefficients: 



a(y,z) 
0(y,z) 

${y,z) 



S23S13 



2B(si3, S23, S123) + ^12(S13, «23, S123) - ^4ll(si3, S23, S123) 



— ^2s 2 3-B(si3, S23, S123) + 2(S12 + Si 3 )Aii(si 3 , S 2 3, S123) 
+ 523(^12(513, «23 ; S123) + ^13(Sl3, S23, Si 23 ))^ , 

^411(513, S23, S123) - ^4l3(si3, S23, S123) 
^4ll(si3, S23, S123)- 



513^23 
4 

S12S13 



(3.8) 



When the hadron tensor is contracted with or the lepton current V 11 , the final term of Eq. (3.7) vanishes^]. 
Furthermore, current conservation implies the following relation between the four helicity coefficients, 



a(y,z) - f3(y,z) -j(y,z) 



2.s 



12:! 



S12 



6(y,z) = 0. 



(3.9) 



This relation is fulfilled automatically once the tensor coefficients are inserted and does therefore not yield a 
further reduction of the tensor basis. 

As with the tensor coefficients, the helicity amplitude coefficients a, (3 and 7 are vectors in colour space 
and have perturbative expansions: 



s And for this reason was omitted in Rcf. 



" (0) + (£)^ ) + (S) 2 ^ )+ ^ 



(3.10) 
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for Q — a, (3, 7. The dependence on (y,z) is again implicit. 

The ultraviolet and infrared properties of the helicity coefficients match those of the tensor coefficients, 











= 57 W 


,un Aq(0),uii 

2e 


fi< 2 ) 







- - 3) Q > (3 - n) 

and 

o« = iW(e)si<® + n( 1 )- finite , 

fi< 2 ) = ^-If«( e )I«( e ) - ^J«( e ) + e- ^~_ 2 ^j + K^j I«(2e) + ff< a >( e )j 0^ 

+I (1) (e)ftW +ft( 2 >< finito , (3.12) 

where J (1) (e) and JT (2) (e) are defined in Eqs. fl2.42j ) and ( |2.44|) respectively. 
At leading order 

a(%,zH/?(°>(y,z) = l and 7 ( %,*) = 0. (3.13) 

The renormalized next-to-leading order helicity amplitude coefficients can be straightforwardly obtained to 
all orders in e from the tensor coefficients using Eqs. ( 2.21 )— ( 2.24 ). For practical purposes, they are needed 
through to C(e 2 ) in evaluating the infrared-divergent one-loop contribution to the two-loop amplitude, while 
only the finite piece is needed for the one-loop self-interference. They can be decomposed according to their 
colour structure as follows: 

n (l), finite ^ z) = Nan{yj z) + }_ 5o(2A z) + fa Cn{yj z) . (3.14) 

The expansion of the coefficients through to e 2 yields HPLs and 2dHPLs up to weight 4 for an, &n and up 
to weight 3 for cq. The explicit expressions are of considerable size, such that we only quote the e°-terms 
here (although these have been known already for a long time The expressions through to C(e 2 ) can be 
obtained in FORM format from the authors. An example of the size and structure of those coefficients can be 
found in ]l8| , where we explicitly list the helicity-averaged one- loop times one-loop and tree times two- loop 
matrix elements. The one-loop coefficients read: 

a a (y, z) = ~-^ + ^H(0; z) - Ih(0, *)G(0; y) - ^H(1,0; z) - ^G(0; y) + ±G(1,0; y) 
' r H(0;z)- ./ . (l + 2H(0;z)) + 0(e) , 



4(1 -z) 2 v ' ' 4(1 -z) 

2 

b a (y,z) = ^(H(0;z)G(l-z;y)+H(l;z)G(-z;y)-G(-z,l-z;y)) +^(-H(0;z) 

+2H(0; z)G(l - z; y) - H(l; z) + 2H(1; z)G(-z; y) + G(l - z;y) - 2G(-z, 1 - z; yjj 

' -H(0; z) - ^-H(0; z) + — 1 H(0; z) + 1 (l + 2H(0; «)) + - - -H(0; z) 



2y(l - z) y ' ' "^y 4(1 - z) 2 v ' ' 4(1 - z) V v ' V 4 4 

+ iH(0; z)G(l -z;y) + ^H(0, 1; z) - jH(l; z) + H(l; z)G(-z; y) - 1h(1; z)G(0; y) 

+ |G(1 - z; y) + ~G(1 - z, 0; y) - G(-z, 1 - z; y) + ~G(0, 1 - z; y) - ^0(1, 0; y) + 0(e) 

1 I 7 7T 

c a (y,z) = --R(0;z)--G(0;y) + -+O(e) , 

*) = -f - ^ + ^H(0; z) - |H(0; z)G(0; y) - ±H(1,0; z) - ^G(0; y) + ±G(1,0; y) 
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4(1 -z) 

My, z) = Z{1 2y2 Z) ( - H(0; z)G(l -z;y)- H(l; z)G(-z; y) + G(-z, 1 - y)) + |- ( - H(0; z) 
+2H(0; z)G(l - z; y) - H(l; z) + 2H(1; z)G(-z; y) + G(l - z; y) - 2G(-z, 1 - z;y) 
+ i- ( - 2H(0; z)G(l -z;y) + H(l; z) - 2H(1; z)G(-z; y) - G(l - z; y) + 2G(-z, 1 - z; y) 

( - H < 15 *> + G(1 - z; y) ) + ^) - W-T) H(0; z) + 2foT3 ( H(1; z) 

-G(l - z; y)) + | - ^H(0; z) + ^(0; z)G(l - z; y) + ^H(0, 1; z) - |h(1; z) 

+H(1; z)G(-z; y) - ^H(l; z)G(0; y) + |g(1 - z; y) + ^G(l - z, 0; y) - G(-z, 1 - z; y) 



+ iG(0,l-z;y)-iG(l,0;y))+O(e), 



I I 7 7T 

C/3 (y,z) = -_H(0;z)--G(0;y) + - + 0(e) , 

<h(y,z) = -i + i ^^H(0;z) + ^ T ^(l-H(0;z))+O( e ) ! 

6 7 (y,z) = \ + ^(- H(0; z)G(l - z; y) - H(l; z)G(-z; y) + G(-z, 1 - z; y)) - - H(0; z) 

+ i ( H (0; z) + H(l; z) - G(l - z; y)) + ( H (l; z) - G(l - z; y)) - 

+ 2(y^) ( - H(1; 2) + G(1 - Z; y) ) - 4(T 1 ^ H(0; Z) + W=T) ( - 1 + H( ° ; Z) ) + ° {6) ' 
c 7 (y,z) = 0. (3.15) 

It should be noted that these finite pieces of the one-loop coefficients can equally well be written in terms of 
ordinary logarithms and dilogarithms, see |],[6j. The reason to express them in terms of HPLs and 2dHPLs 
here is their usage in the infrared counter-term of the two-loop coefficients, which cannot be fully expressed 
in terms of logarithmic and poly logarithmic functions. 

The finite two-loop remainder is obtained by subtracting the predicted infrared structure (expanded 
through to C(e )) from the renormalized helicity coefficient. We further decompose the finite remainder 
according to the colour structure, as follows: 

ft( 2 >> finitc (y,z) = N 2 An(y,z) + Bn(y,z) + ^Cn(y,z) + NN F D n (y,z) 

+^E n (y, z) + N 2 F F n (y, z) + N F , V (A - n) G (y, z) , (3.16) 

where the last term is generated by graphs where the virtual gauge boson does not couple directly to the 
final-state quarks. This contribution is denoted by N F ,v and is proportional to the charge weighted sum of 
the quark flavours. In the case of purely electromagnetic interactions we find, 

N F ^ = ^ifi . (3.17) 

Including Z-interactions, the same class of diagrams yields not only a contribution from the vector component 
of the Z, which for the right-handed quark amplitude is given by 

Nf.z = E? (L 2 | Z + ^ , (3.18) 
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but also a contribution involving the axial couplings of the Z |64[] . This contribution vanishes if summed over 
isospin doublets. The large mass splitting of the third quark family induces a non-vanishing contribution 
from this class of diagrams, which can however not be computed within the framework of massless QCD 
employed here, but can only be obtained within an effective theory with large top-quark mass. In contrast to 
the vector contribution from these diagrams, which is finite, one could expect divergences in the axial vector 
contribution, which would be cancelled by the single unresolved limits of the corresponding axial contributions 
to four-parton final states [p^ , |l5| . Results from the four-parton final states show that this axial contribution 
is numerically very small |65| . 

The helicity coefficients contain HPLs and 2dHPLs up to weight 4 in the A, B, C, G-terms, up to weight 3 
in the D, i?-terms (which do moreover contain only a limited subset of purely planar master integrals) and up 
to weight 2 in the F-term. The size of each helicity coefficient is comparable to the size of the helicity-averaged 
tree times two-loop matrix element quoted in ||l8| . We do therefore only quote the A- and Z?-terms of each 
coefficient, which form the leading colour contributions, and which turn out to be numerically dominant, 
approximating the full expressions to an accuracy of about 20%. The complete set of coefficients in FORM 
format can be obtained from the authors. 

These leading colour terms are: 



A a {y,z) = 
1 



48y(l 
+6G(l,0;y) 

+5G(l,0;y) 

-28G(l,0;y) 
1 



7T 2 - 13H(0; z) + 6H(1, 0; z) + 6G(1, 0; y) 



48|/(l-y-z) 



vr 2 -13H(0;z) + 6H(l,0;z) 



16(1 - y y 

z 



:G(0;y) 



16(1 -y) 12(1 -y-zf 



~6~ 



16(1 -y- z) 



Utt 2 



+ 



G(0;y) 



16(1 -y-zf 
1 



[1-zf 



G 



5H(0;z)G(0 ;2 ,)-5H(l,0;z) 
- 11H(0; z) + 28H(0; z)G(0; y) + 28H(1, 0; z) + 11G(0; y) 
■ + 11H(0; z)G(0; y) + 11H(1, 0; z) - 11G(1, 0; y) 
+ ir 2 (3H(0; z) + 3H(1; z) - G(l -z;y) + G(0; y)) + 6C 3 



3(1 -y) 

^-H(0; z) - 6H(0; z)G(l - z, 0; y) + 10H(0; z)G(0; y) + 45H(0, 0; z) + 12H(0, 0; z)G(0; y) 

+18H(0, 1,0; z) - H(l, 0; z) - 6H(1, 0; z)G(l -z;y) + 6H(1, 0; z)G(0; y) + 12H(1, 0, 0; z) + 18H(1, 1,0; z) 

+6G(l-z,l,0;j/)-6G(0,l,0;y) ' 
277 



72(1 - z) 



7T 2 ( - 8 + 9H(0; z) + 9H(1; z) - 3G(1 - z; y) 



+3G(0; y)) + I8C3 - ^ - 65H(0; z) - 18H(0; z)G(l - z, 0; y) + 39H(0; z)G(0; y) + 81H(0, 0; z) 
+36H(0, 0; z)G(0; y) + 54H(0, 1, 0; z) - 48H(1, 0; z) - 18H(1, 0; z)G(l -z;y) + 18H(1, 0; z)G(0; y) 
+36H(1, 0, 0; z) + 54H(1, 1, 0; z) + 18G(1 - z, 1, 0; y) + 15G(0; y) - 18G(0, 1, 0; y) - 9G(1, 0; y) 

' 2tt 2 H(1; z) - tt 2 G(0; y) + 12C 3 - 6H(1, 0; z)G(0; y) - 12H(1, 1, 0; z) - 6G(0, 1, 0; y) 



48(1 -y-zf 
1 



48(1 -y-z) 

-20G(0;y)-12G(0,l,0;y) 



4tt 2 H(1; z) - 2tt 2 G(0; y) + 24C 3 - 13H(0; z) - 12H(1, 0; z)G(0; y) - 24H(1, 1, 0; z) 
2 - 928 



7T 

288 



— - 5H(0; z) + 12H(0; z)G(l -z;y) + 36H(0; z)G(0; y) 



-12H(0; z)G(l; y) + 24H(0, 1; z) + 24H(1; z)G(l - z\y) - 24H(1; z)G(-z; y) - 12H(1; z)G(l; y) 
-24H(1, 0; z) + 12H(1, 1; z) - 44G(1 - z; y) + 12G(1 - z, 0; y) - 24G(1 - z, l;y) + 24G(-z, 1 - 2; y) 

-24G(0, 1 - z; y) + 49G(0; y) - 24G(0, l;y)+ 12G(1, 1 - z; y) + 44G(1; y) - 36G(1, 0; y) + 24G(1, 1; y) 



72 



317 - 18H(0; z) + 90H(1; z) - 72G(1 - z; y) - 18G(0; y) - 18G(1; y) 



IItt 4 
360 



1 

72 



89959 
144 



14 



2149 
12 



H(0; z) - 66H(0; z)G(l - z, 0; y) - 18H(0; z)G(l - z, 1, 0; y) + 36H(0; z)G(-z, l-z,0; y) 



-36H(0; z)G(0, l-z,0;y)- 66H(0; z)G(0; y) + 126H(0; z)G(0, 0; y) - 18H(0; z)G(0, 1, 0; y) 

23 

+18H(0; z)G(l, 1 - z, 0; y) - 3H(0; z)G(l, 0; y) - 36H(0; z)G(l, 0, 0; y) + yH(0, 0; z) 

+72H(0, 0; z)G(0; y) + 36H(0, 0; z)G(0, 0; y) + 72H(0, 0, 1, 0; z) + 3H(0, 1, 0; z) - 18H(0, 1, 0; z)G(l - z; y) 

+36H(0, 1, 0; z)G(-z; y) + 18H(0, 1, 0; z)G(0; y) - 18H(0, 1, 0; z)G(l; y) + 36H(0, 1, 1, 0; z) 

-71H(1, 0; z) - 66H(1, 0; z)G(l - z; y) + 18H(1, 0; z)G(l - z, 0; y) + 36H(1, 0; z)G(-z, 1 - z; y) 

-36H(1, 0; z)G(-z, 0; y) - 36H(1, 0; z)G(0, 1 - z; y) + 96H(1, 0; z)G(0; y) + 18H(1, 0; z)G(l, 1 - z; y) 

-18H(1, 0; z)G(l, 0; y) + 72H(1, 0, 0; z) + 36H(1, 0, 0; z)G(0; y) + 72H(1, 0, 1, 0; z) 

+36H(1, 1, 0; z)G(l - z; y) - 36H(1, 1, 0; z)G(-z; y) - 18H(1, 1, 0; z)G(l; y) + 36H(1, 1, 0, 0; z) 

+18H(1, 1, 1,0; z) + 18G(1 - z,0, 1,0; y) + 66G(1 - z, 1,0; y) + 36G(1 - z, 1, 1,0; y) 

49 

-36G(-z, 1 - z, 1, 0; y) - 36G(-z, 0, 1, 0; y) + 36G(0, 1 - z, 1, 0; y) + — G(0; y) + 160G(0, 0; y) 
-36G(0, 0, 1, 0; y) - 30G(0, 1, 0; y) + 36G(0, 1, 1, 0; y) - 18G(1, 1 - z, 1, 0; y) + 71G(1, 0; y) 
-126G(1, 0, 0; y) + 54G(1, 0, 1, 0; y) - 66G(1, 1, 0; y) + 36G(1, 1, 0, 0; y) - 36G(1, 1, 1, 0; y) 



11 



r H(0;z) 



11 - 22H(0;z) 



16(1 - z) 2 " v "'~ y ' 16(1 -z) 
-66H(0; z)G(0; y) - 66H(1, 0; z) - 110G(0; y) + 66G(1, 0; y) 



44tt 2 2345 



18 



llH(0;z) 



D a {y,z) 



12y(l-z) 
-G(l,0;y) 

z 2 

+ 



H(0;z) 



12y(l-y-z) 



H(0;z) + 



6(1 -y-zf 



+ H(0;z)G(0;y)+H(l,0;z) 



+ 



4(1 -y-z) 



— + H(0; z) - 2H(0; z)G(0; y) - 2H(1, 0; z) - G(0; y) + 2G(1, 0; y) 



4(1 -y-zf 
1 



— - H(0; z)G(0; y) - H(l, 0; z) + G(l, 0; y) 



1 



12(1 -y) 



G(0;y) 



+ 



72(1 -z) 2 
1 



7T 2 + 25H(0; z) - -H(0; z)G(0; y) - 9H(0, 0; z) + 6H(1, 0; z) 



144(1 - z) 
1 



4tt 2 + 38 + 37H(0; z) - 6H(0; z)G(0; y) - 36H(0, 0; z) + 24H(1, 0; z) - 3G(0; y) 



12(1 -y-z) 

19, 1 

C3 H 

36^ 144 



H(0;z) + 2G(0;y) 



72 



395 

— - H(0; z) + 2G(1 - z; y) - G(0; y) - 2G(1; y) 



3661 



18 



25H(0; z) + 24H(0; z)G(l - z, 0; y) + 29H(0; z)G(0; y) - 36H(0; 2)G(0, 0; y) 



+6H(0; z)G(l, 0; y) - 28H(0, 0; z) - 36H(0, 0; z)G(0; y) - 6H(0, 1, 0; z) + 40H(1, 0; z) 

+24H(1, 0; z)G(l - z; y) - 30H(1, 0; z)G(0; y) - 36H(1, 0, 0; z) - 24G(1 - z, 1, 0; y) + 53G(0; y) 

-82G(0, 0; y) + 6G(0, 1, 0; y) - 40G(1, 0; y) + 36G(1, 0, 0; y) + 24G(1, 1, 0; y) 



+Z7T< 



r H(0;z) 



1 + 2H(0; z) 



8(1 -z) 2 "^'"' ' 8(1 -z) 
+ 12H(l,0;z) + 31G(0;y)- 12G(l,0;y) 



1 

48 



8tt 2 28 



13H(0;z) + 12H(0;z)G(0;y) 



Af){y,< 



15 



16(1 -y) 2 



G(0;y) 



+ 



16(1 -y) 16(y + z) 2 



^-H(l; z) - ^-G(l -*;»)- 94H(0; z)G(l - z, 0; y) 



-94H(0, 1, 0; z) - 99H(1, 0; z) - 94H(1,0; z)G(l - z;y) + 94H(1, 0; z)G(0; y) + 94H(1, 1, 0; z) 



+94G(1 - z, 1, 0; y) + 94G(0, 1, 0; y) - 99G(1, 0; y) 



+ 



A7tt 2 



16(y + z)L 

-94H(0; z)G(0; y) - 94H(1, 0; z) - 55G(0; y) + 94G(1, 0; y) 



12(1 -y-z) 2 



+ 11 + 44H(0;z) 



5 " 2 71-2 HC1 • 7 s ! 



~G(0; y) + 3C 3 + 5H(0; z)G(0; y) + 5H(1, 0; z) - |h(1, 0; z)G(0; y) - 3H(1, 1, 0; z) - |g(0, 1, 0; y) 



-5G(l,0;y) 



+38G(l,0;y) 



12(1 -y-z) 

^2 



197T 2 53 
- — + 5H(0; z) - 38H(0; z)G(0; y) - 38H(1, 0; z) - T G(0; y) 

11tt 2 H(1; z) + 11tt 2 G(1 -z;y) + 66H(0; z)G(l - z, 0; y) + 66H(0, 1,0; z) 



8(y + z)3 

+33H(1, 0; z) + 66H(1, 0; z)G(l - z;y) - 66H(1, 0; z)G(0; y) - 66H(1, 1, 0; z) - 66G(1 - z, 1, 0; y) 



-66G(0,l,0;y) + 33G(l,0;y) 
+33G(0;y)-132G(l,0;y) 
+HG(0;y)-66G(l,0;y) 



16(y + z) 2 



+ 22tt 2 - 33H(0; z) + 132H(0; z)G(0; y) + 132H(1, 0; z) 



16(y + z) 

~2 



IItt 2 - 11H(0; z) + 66H(0; z)G(0; y) + 66H(1, 0; z) 



+ 



48(1 -y-z) 2 



53tt 2 



- 53H(0; z)G(0; y) - 53H(1, 0; z) + 53G(1, 0; y) 



+ 



16(1 — y — z 

.3 rll7r 2 



IItt 2 - 11H(0; z) + 66H(0; z)G(0; y) + 66H(1, 0; z) + 11G(0; y) - 66G(1, 0; y) 



IItt 2 

H(l; z) — G(l - z; y) - 33H(0; z)G(l - z, 0; y) - 33H(0, 1, 0; z) 



8(y + z) 4 

-33H(1, 0; z)G(l - z; y) + 33H(1, 0; z)G(0; y) + 33H(1, 1, 0; z) + 33G(1 - z, 1, 0; y) + 33G(0, 1, 0; y) 



+ 



8(y + z)3 



ll^ 2 



33H(0; z)G(0; y) - 33H(1, 0; z) + 33G(1, 0; y) 



16(y + z) 2 



IItt 2 



-33H(0; z)G(0; y) - 33H(1, 0; z) + 33G(1, 0; y) 
.3 rll7r 2 



8(y + z) 



IItt 2 



- 11H(0; z)G(0; y) - 11H(1, 0; z) 



+HG(l,0;y) 

,3 



+ 



8(1 -y-z) 
1 



16(1 -y-z) 2 
ll^ 2 



+ 11H(0; z)G(0; y) + 11H(1, 0; z) - 11G(1, 0; y) 



6 



11H(0; z)G(0; y) - 11H(1, 0; z) + 11G(1, 0; y) 



+ 



24(1 -y) 



G(0;y) 



i 

tt 2 (- - 3H(0; z) - 3H(1; z) + G(l - z; y) - G(0; y)) - 6Cs + — H(0; z) 



48(1 - z) 

+6H(0; z)G(l - z, 0; y) - 10H(0; z)G(0; y) - 45H(0, 0; z) - 12H(0, 0; z)G(0; y) - 18H(0, 1, 0; z) 
+H(1, 0; z) + 6H(1, 0; z)G(l - z; y) - 6H(1, 0; z)G(0; y) - 12H(1, 0, 0; z) - 18H(1, 1, 0; z) 

1 



-6G(l-z,l,0;y) + 6G(0,l,0;y) 



+ 



^H(l; z) + ^G(l - z; y) + 14H(0; z)G(l - z, 0; y) 



8(y + ^) 

+14H(0, 1, 0; z) + 25H(1, 0; z) + 14H(1, 0; z)G(l - z; y) - 14H(1, 0; z)G(0; y) - 14H(1, 1, 0; z) 

X r^2 ^2 



-14G(1 - z, 1, 0; y) - 14G(0, 1, 0; y) + 25G(1, 0; y) 



8(1 -y-«) s 



2 

r ^-H(l;z) + ^G(0;y)-2C 3 



+H(1, 0; z)G(0; y) + 2H(1, 1, 0; z) + G(0, 1, 0; y) 



+ 



24(1 -y-z) 



13^ 2 



+ 10H(0;z)G(0;y) + 13H(l,0;z) 



16 



+10G(0;y)-7G(l,0;y) 



7T 

288 



2 r 04g 

- — - 5H(0; z) + 12H(0; z)G(l -z;y) + 36H(0; z)G(0; y) 



-12H(0; z)G(l; y) + 24H(0, 1; z) + 24H(1; z)G(l - z; y) - 24H(1; z)G(-z; y) - 12H(1; z)G(l; y) 
+24H(1, 0; z) + 12H(1, 1; z) - 44G(1 -z;y) + 12G(1 - z, 0; y) - 24G(1 - z, 1; y) + 24G(-z, 1 - z; y) 

-24G(0, 1 - z; y) + 49G(0; y) - 24G(0, 1; y) + 12G(1, 1 - z; y) + 44G(1; y) - 36G(1, 0; y) + 24G(1, 1; y) 



IItt 4 Cs 

+ ^7T + 



+ 



1 



144 



| _ 2 L 317 - 18H(0; z) + 90H(1; z) - 72G(1 - z; y) - 18G(0; y) - 18G(1; y) 
1735 

H(0; z) - 132H(0; z)G(l - z, 0; y) - 36H(0: z)G(l - z, 1, 0: y) + 72H(0; z)G(-z, 1 - z, 0; y) 

6 



79987 



72 



-72H(0; z)G(0, 1 - z, 0; y) - 150H(0; z)G(0; y) + 252H(0; z)G(0, 0; y) - 36H(0; z)G(0, 1, 0; y) 

+36H(0; z)G(l, 1 - z, 0; y) - 6H(0; z)G(l, 0; y) - 72H(0; z)G(l, 0, 0; y) + 23H(0, 0; z) 

+144H(0, 0; z)G(0; y) + 72H(0, 0; z)G(0, 0; y) + 144H(0, 0, 1, 0; z) + 6H(0, 1, 0; z) 

-36H(0, 1, 0; z)G(l - z; y) + 72H(0, 1, 0; z)G(-z; y) + 36H(0, 1, 0; z)G(0; y) - 36H(0, 1, 0; z)G(l; y) 

+72H(0, 1, 1, 0; z) - 160H(1, 0; z) - 132H(1, 0; z)G(l - z; y) + 36H(1, 0; z)G(l - z, 0; y) 

+72H(1, 0; z)G(-z, 1 - z; y) - 72H(1, 0; z)G(-z, 0; y) - 72H(1, 0; z)G(0, 1 - z; y) + 192H(1, 0; z)G(0; y) 

+36H(1, 0; z)G(l, 1 - z; y) - 36H(1, 0; z)G(l, 0; y) + 144H(1, 0, 0; z) + 72H(1, 0, 0; z)G(0; y) 

+144H(1, 0, 1, 0; z) + 72H(1, 1, 0; z)G(l - z; y) - 72H(1, 1, 0; z)G(-z; y) - 36H(1, 1, 0; z)G(l; y) 

+72H(1, 1,0, 0; z) + 36H(1, 1, 1, 0; z) + 36G(1 - z, 0, 1, 0; y) + 132G(1 - z, 1, 0; y) + 72G(1 - z, 1, 1, 0; y) 

-72G(-z, 1 - z, 1, 0; y) - 72G(-z, 0, 1, 0; y) + 72G(0, 1 - z, 1, 0; y) + ^G(0; y) + 320G(0, 0; y) 
-72G(0, 0, 1, 0; y) - 60G(0, 1, 0; y) + 72G(0, 1, 1, 0; y) - 36G(1, 1 - z, 1, 0; y) + 160G(1, 0; y) 
-252G(1, 0, 0; y) + 108G(1, 0, 1, 0; y) - 132G(1, 1, 0; y) + 72G(1, 1, 0, 0; y) - 72G(1, 1, 1, 0; y) 



+ Hw^) H(0; ' )+2C3+ ^ 



44tt 2 1751 
~~3 18~ 



11H(0; z) - 66H(0; z)G(0; y) - 66H(1, 0; z) 



-110G(0;y) + 66G(l,0;y) 



Dp(y,z) = 

z 



4tt 2 4tt 2 

— H(l; z) + — G(l - z; y) + 8H(0; z)G(l - z, 0; y) + 8H(0, 1, 0; z) + 9H(1, 0; z) 



4(y + z) 2 

+8H(1, 0; z)G(l - z; y) - 8H(1, 0; z)G(0; y) - 8H(1, 1, 0; z) - 8G(1 - z, 1, 0; y) - 8G(0, 1, 0; y) 



+9G(l,0;y) 

z 

+ 



4(y + z) 



^- - 1 - 4H(0; z) + 8H(0; z)G(0; y) + 8H(1, 0; z) + 5G(0; y) - 8G(1, 0; y) 



6(1 -y-z) 2 



H(0;z)G(0;y)-H(l,0;z) + G(l,0;y) 



+ 



12(1 -y-z) 



7^ 2 



2H(0;z) 



+14H(0; z)G(0; y) + 14H(1, 0; z) + 5G(0; y) - 14G(1, 0; y) 



tt 2 H(1;z) -7T 2 G(l-z;y) 



2(y + z)3 

-6H(0; z)G(l - z, 0; y) - 6H(0, 1,0; z) - 3H(1, 0; z) - 6H(1, 0; z)G(l -z;y) + 6H(1, 0; z)G(0; y) 

^2 



+6H(1, 1, 0; z) + 6G(1 - z, 1, 0; y) + 6G(0, 1, 0; y) - 3G(1, 0; y) 
-12H(0; z)G(0; y) - 12H(1, 0; z) - 3G(0; y) + 12G(1, 0; y) 



+ 



4(y + z) 2 



2tt 2 + 3H(0;z) 



-6H(0; z)G(0; y) - 6H(1, 0; z) - G(0; y) + 6G(1, 0; y) 



4(y + z) 

^2 



+ 



12(1 -y-z) 2 



-^ 2 + H(0;z) 
5tt 2 



+ 5H(0;z)G(0;y) 



17 



4(1 -y-z) 



7T 2 + H(0; z) - 6H(0; z)G(0; y) - 6H(1, 0; z) - G(0; y) 



+5H(l,0;z)-5G(l,0;y) r 
+3H(1, 0; z)G(l - «; y) - 3H(1, 0; z)G(0; y) - 3H(1, 1, 0; z) - 3G(1 - z, 1, 0; y) - 3G(0, 1, 0; y) 



2 2 

yH(l; z) + yG(l -z;y) + 3H(0; z)G(l - z, 0; y) + 3H(0, 1, 0; z) 



2(y + z) 3 
+3H(l,0;z)-3G(l,0;y) 



y + 3H(0; z)G(0; y) + 3H(1, 0; z) - 3G(1, 0; y) 



4(y + z) 



y + 3H(0;z)G(0;y) 



2(y + z) 



•H(0;«)G(0;y)+H(l,0;«)-G(l,0;y) 



+ 



- — - H(0; z)G(0; y) - H(l, 0; z) + G(l, 0; y) 



2(1 -y-z) 



-+H(0;z)G(0;y) 



1 



12(1 -y) 
1 



G(0;y) 



72(1 - z) 



^ 2 -25H(0;z) + -H(0;z)G(0;y) 



4(1 -y-z) 2 
+H(l,0;z)-G(l,0;y) 

+9H(0,0;z)-6H(l,0;z) , 

J 2(y + z) 

-2H(1, 0; z) - H(l, 0; z)G(l - z; y) + H(l, 0; z)G(0; y) + H(l, 1, 0; z) + G(l - z, 1, 0; y) + G(0, 1, 0; y) 



2 2 

yH(l; z) - 1-G(1 - z; y) - H(0; z)G(l - z, 0; y) - H(0, 1, 0; z) 



-2G(l,0;y) 
- 2 r395 



6(1 -y-z) 



H(0; z)G(0; y) - H(l, 0; z) - G(0; y) + G(l, 0; y) 



7T 

72 



12 



- H(0; z) + 2G(1 - z; y) - G(0; y) - 2G(1; y) 



19, 1 

— C3 H 

36 s 144 



2977 



18 



- 10H(0;z) 



+24H(0; z)G(l - z, 0; y) + 29H(0; z)G(0; y) - 36H(0; z)G(0, 0; y) + 6H(0; z)G(l, 0; y) - 28H(0, 0; z) 
-36H(0, 0; z)G(0; y) - 6H(0, 1, 0; z) + 40H(1, 0; z) + 24H(1, 0; z)G(l - z; y) - 30H(1, 0; z)G(0; y) 
-36H(1, 0, 0; z) - 24G(1 - z, 1, 0; y) + 56G(0; y) - 82G(0, 0; y) + 6G(0, 1, 0; y) - 40G(1, 0; y) 

+36G(l,0,0;y) + 24G(l,l,0;y) 

+ H ~ 8 (l-z) H(0; Z) + i " f + 13H( ° ; Z) + 12H( ° ; Z)G( °' y) + 12H(1 ' ° 5 Z) + 31G( ° ; y) 



-12G(l,0;y) 



A 7 (y,z) 



1 



48y(l - z) 
+6H(l,0;z) + 6G(l,0;y) 



7T 2 + 13H(0; z) - 6H(1, 0; z) - 6G(1, 0; y) 



16(y + zf 



lln 2 



48y(l -y-z) 
ll^ 2 



7T 2 - 13H(0;«) 



H(l; z) + -— G(l - z; y) + 22H(0; z)G(l - z, 0; y) 



+22H(0, 1,0; z) + 55H(1, 0; z) + 22H(1, 0; z)G(l - z; y) - 22H(1, 0; z)G(0; y) - 22H(1, 1, 0; z) 



-22G(1 - z, 1, 0; y) - 22G(0, 1, 0; y) + 55G(1, 0; y) 



+ 



16(y + z) 

+22H(0; z)G(0; y) + 22H(1, 0; z) + 33G(0; y) - 22G(1, 0; y) 



IItt 2 



-2C 3 + H(l, 0; z)G(0; y) + 2H(1, 1, 0; z) + G(0, 1, 0; y) 



+ 



8(1 -y-zY 
10^ 2 



11 - 22H(0;z) 



yH(l;z) + ^G(0;y) 



48(1 -y-z) 



+20H(0; z)G(0; y) + 20H(1, 0; z) + 20G(0; y) - 20G(1, 0; y) 
22tt 2 



8(y + z)3 



3 

22tt 2 



+ 13H(0;z) 
H(l;z) 



-G(l - z; y) - 44H(0; z)G(l - z, 0; y) - 44H(0, 1, 0; z) - 33H(1, 0; z) - 44H(1, 0; z)G(l - z; y) 



18 



+44H(1, 0; z)G(0; y) + 44H(1, 1, 0; z) + 44G(1 - z, 1, 0; y) + 44G(0, 1, 0; y) - 33G(1, 0; y) 



16(y + z) ; 

,2 



16(y + z) 

,2 



44tt 2 
~ ~~ 3~ 
22n 2 
3 

57T 



12(1 - 2/ - z) 5 



+ 33H(0; z) - 88H(0; z)G(0; y) - 88H(1, 0; z) - 33G(0; y) + 88G(1, 0; y) 
- 11H(0; z) - 44H(0; z)G(0; y) - 44H(1, 0; z) - 11G(0; y) + 44G(1, 0; y) 

+ 



— + 5H(0; z)G(0; y) + 5H(1, 0; z) - 5G(1, 0; y) 
o 



16(1 -y-z) 



22^ 2 



8(y + z) 4 



+11H(0; z) - 44H(0; z)G(0; y) - 44H(1, 0; z) - 11G(0; y) + 44G(1, 0; y)J + 
IItt 2 

+— G(l - z; y) + 33H(0; z)G(l - z, 0; y) + 33H(0, 1, 0; z) + 33H(1, 0; z)G(l - z; y) 



-33H(1, 0; z)G(0; y) - 33H(1, 1, 0; z) - 33G(1 - z, 1, 0; y) - 33G(0, 1, 0; y) 

.3 rU?r 2 



+ 



8(y + z)3 



3 

llTT 2 



llTT 2 



H(l;z) 



+33H(0; z)G(0; y) + 33H(1, 0; z) - 33G(1, 0; y) 



16(y + z) 2 



+ 33H(0; z)G(0; y) + 33H(1, 0; z) 



+ 



-33G(l,0;y) 

z 3 

~16(1- y-z) 2 



.3 rll7r 2 



8(y + z) 

" llTT 2 



11H(0; z)G(0; y) + 11H(1, 0; z) - 11G(1, 0; y) 



11H(0; z)G(0; y) - 11H(1, 0; z) + 11G(1, 0; y) 



+ 



8(1 -y-z) 



IItt 2 



+ llH(0;z)G(0;y) + HH(l,0;z) - HG(l,0;y) 
355, 



+ 



48(1 - z) 2 



^ 2 (--3H(0;z) 



-3H(1; z) + G(l - z; y) - G(0; yj) - 6(3 + — H(0; z) + 6H(0; z)G(l - z, 0; y) - 10H(0; z)G(0; y) 
-45H(0, 0; z) - 12H(0, 0; z)G(0; y) - 18H(0, 1, 0; z) + H(l, 0; z) + 6H(1, 0; z)G(l - z; y) 
-6H(1, 0; z)G(0; y) - 12H(1, 0, 0; z) - 18H(1, 1, 0; z) - 6G(1 - z, 1, 0; y) + 6G(0, 1, 0; y) 



+ 



1 



48(1 - z) 



7 



277 493. 



tt 2 ( - - + 3H(0; z) + 3H(1; z) - G(l - z; y) + G(0; y)) + 6(3 + — — H(0; z) 



-6H(0; z)G(l - z, 0; y) + 4H(0; z)G(0; y) + 45H(0, 0; z) + 12H(0, 0; z)G(0; y) + 18H(0, 1, 0; z) 
-7H(1, 0; z) - 6H(1, 0; z)G(l -z;y) + 6H(1, 0; z)G(0; y) + 12H(1, 0, 0; z) + 18H(1, 1, 0; z) 

1 



+6G(1 - z, 1, 0; y) - 10G(0; y) - 6G(0, 1, 0; y) + 6G(1, 0; y) 



+ 



48(1 -y-z) 



- 7T 2 + 13H(0; z) 



-6H(l,0;z)-6G(l,0;y) 



1 

48 
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+ 23H(0; z) + 6H(0; z)G(0; y) + 6H(1, 0; z) + 10G(0; y) 



-6G(l,0;y) 



+ ITT< 



lL_ H (0;z) + ^l 



16(1 - z) 



16(1 - z) 



1-H(0;z) 



11 
16 



D 7 (y;z) = 



1 



-H(0;z) 



1 



-H(0;z) 



_H(l;z)- y G(l-z;y) 



12y(l - z)" v "'~ y ' 12y(l-y-z)" v ~'~ y ' 4(y + z) 2 

-2H(0; z)G(l - z, 0; y) - 2H(0, 1,0; z) - 5H(1, 0; z) - 2H(1, 0; z)G(l - z; y) + 2H(1, 0; z)G(0; y) 

^2 



+2H(1, 1,0; z) + 2G(1 - z, 1,0; y) + 2G(0, 1,0; y) - 5G(1, 0; y) 
-2H(0; z)G(0; y) - 2H(1, 0; z) - 3G(0; y) + 2G(1, 0; y) 



4(y + z) 



-2H(0; z)G(0; y) - 2H(1, 0; z) - 2G(0; y) + 2G(1, 0; y) 



+ 



12(1 -y-z) 
z 2 

2(y + z)3 



- — + l + 2H(0;z) 

7T 2 

y-H(0;z) 



2tt 2 TT/1 . 27r 2 ^. . 
— H(l;z) + — G(l -z;y) 
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+4H(0; z)G(l - z, 0; y) + 4H(0, 1, 0; z) + 3H(1, 0; z) + 4H(1, 0; z)G(l -z;y)- 4H(1, 0; z)G(0; y) 

i z 2 r47r 2 

-4H(1, 1, 0; *) - 4G(1 - *, 1, 0; y) - 4G(0, 1, 0; y) + 3G(1, 0; y) 



+8H(0; z)G(0; y) + 8H(1, 0; z) + 3G(0; y) - 8G(1, 0; y) 

-4H(l,0;z) + G(0;y)-4G(l,0;y)' 
- 2 -2ir 2 



4(y + z) 



4(y + z) 2 
3 



— -3H(0;z) 
H(0;z) + 4H(0;z)G(0;y) 



6(1 -y-zf 



H(0;z)G(0;y)-H(l,0;z) + G(l,0;y) 



4(1-1/-*) 

■v3 



— - H(0; z) + 4H(0; z)G(0; y) + 4H(1, 0; z) + G(0; y) ~ 4G(1, 0; y) 



2 2 

y H(l; z) - yG(l - z; y) - 3H(0; z)G(l - z, 0; y) - 3H(0, 1, 0; z) - 3H(1, 0; z)G(l - z; y) 



2(y + « 

+3H(1, 0; z)G(0; y) + 3H(1, 1, 0; z) + 3G(1 - z, 1, 0; y) + 3G(0, 1, 0; y) 



2(y + z)3 



Y 



-3H(0; z)G(0; y) - 3H(1, 0; z) + 3G(1, 0; y) 



4(y + z) 2 



— -3H(0;z)G(0;y)-3H(l,0;z) 



-3G(l,0;y) 



2(y + z) 



G 



-H(0;z)G(0;y)-H(l,0;z) + G(l,0;y) 



4(1 -y-z) 2 



7T 

~6~ 



+H(0; z)G(0; y) + H(l, 0; z) - G(l, 0; y) 



2(1 -y-z) 



G 



H(0;z)G(0;y)-H(l,0;z) + G(l,0;y) 



+ 



1 



- 7T 2 - 25H(0; z) + -H(0; z)G(0; y) + 9H(0, 0; z) - 6H(1, 0; z) 



72(1 -z) 2 

65H(0; z) - 3H(0; z)G(0; y) - 18H(0, 0; z) + 12H(1, 0; z) + 3G(0; y) 

1 , 1 



144(1 - z) 



1 



2tt 2 - 38 
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12(1 -y-z) H ^ + 72 



K 1 

— H(0; z) G(0; y) + mt^ ; - Nn H(0; z) + — 

48 v ' ' 48 v ,yj | 8(1 -z) 2 v ' y 8(1 -z) 



1 +H(0;z; 



(3.19) 



From the ftW.nnite and n (2),finite 

it is possible to recover the finite pieces of the helicity-averaged tree 
times two- loop and one-loop squared matrix elements by squaring (3.7): 



Finite (2x0) (x,y,z) 



8VK 



(l-y)(l-y-z) a{2)finitc ^ z) + l^y_ pi 2)Mte {yjZ) 



- 7 (2) ^ nite (y,z) + (y^z) 



J 7 inite {lxl) (x 1 y,z) = 4V K 



(i- y-z) 



(i-y-z) l 

2 



yz 



(1). finite 



1 — y — z z 



+ -) \^ aite (y,z)\ 2 + 
yj 



(y,z)\ 2 

i-y-z + y _\ wl>Al]tii ,. {! ^ ][2 



-3 + y + z 



2- 2z 



2 : 
(1) ' finito (y,z)/3* (1) ' finitc (y,z) 



-(1-y-z) a«> fi " ite (y,*) 7 *«> finite (y,z) 

-(1 + y + z) /3«' finite (y, z) 7 *M' finite (y, z) + (y ^ z) 



(3.20) 



It is important to notice that ( |3.20 ) corresponds, by the very nature of the Weyl-van der Waerden helicity 
formalism, to a scheme with external momenta and polarization vectors in four dimensions (internal states are 
always taken to be d-dimensional), which is sometimes called the 't Hooft-Veltman scheme (28). This scheme 
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is different from the conventional dimensional regularization used in |l8f|, w here all external momenta and 



polarization vectors are (i-dimensional. Nevertheless, one obtains from ( 3.2C ) the same J r inite^ 2x0 \x, y, z) as 
in Eq. (4.17) and Tinite^ lxl \x 1 y, z) as in Eq. (4.25) of fl8fl , since all scheme-dependent terms are correctly 
accounted for by the finite contributions arising from expanding the tree level and one-loop contributions in 
the renormalization and infrared factorization formulae. 

It should also be noted that only the O(e ) terms of J^ 1 ).* 51114 * 3 contribute to !Finite^ lxl \x,y, z), terms 
subleading in e are not required, since no term is multiplied with a divergent factor. Comparing the size 



of these O(e ) terms ( 3.15 ) with the size of J r inite^ lxl \x,y,z) in flq| , it becomes clear that the squared 
one-loop amplitude can be evaluated much more elegantly by squaring the finite remainders of the hclicity 
amplitudes than by computing the squared matrix element. 

4 Conclusions and Outlook 

In this paper, we have presented analytic formulae for the one- and two- loop virtual hclicity amplitudes to the 
process 7* — ► qqg. These amplitudes have been derived by defining projectors, which isolate the coefficients 
of the most general tensorial structure of the matrix element at any order in perturbation theory. Once 
the general tensor is known, the helicity amplitudes follow in a straightforward manner - they are linear 
combinations of the tensor coefficients. We applied the projectors directly to the Feynman diagrams and used 
the conventional approach of relating the ensuing tensor integrals to a basis set of master integrals. This latter 
step is identical to that employed to evaluate the interference of tree- and two-loop graphs in Ref. p8| , apart 
from the fact that the projector is no longer the tree-level amplitude. As anticipated, the finite remainder 
from the interference of tree- and two-loop amplitudes can be reconstructed from the appropriate helicity 
amplitudes, with the difference between treating the external states in d dimensions or four dimensions being 
isolated in the infrared-singular terms. 

The results presented here therefore complement the earlier calculation of the interference of tree- and two- 
loop graphs in Ref. [Q . Knowledge of the helicity amplitudes allows additional information on the scattering 
process. In particular, observables that require knowledge of the polarization tensor of the virtual photon, 
such as oriented event shapes in unpolarized e + e~ scattering or event shapes in polarized e + e~ scattering, 
can be described at two-loop order. 

Similar results can in principle be obtained for (2 + l)-jet production in deep inelastic ep scattering or 
(V + l)-jet production in hadron-hadron collisions. However, the rather different domains of convergence of 
the HPLs and 2dHPLs makes this a non-trivial task, which is discussed in a separate paper |}(|. Nevertheless, 
the helicity approach will provide information on the direction of the decay leptons in (V + l)-jet production 
(with or without polarized protons). Determination of the polarized parton distribution functions in polarized 
electron-proton scattering will also benefit from the knowledge of the two-loop helicity amplitudes in the 
appropriate kinematic region. 

Even though the evaluation of two-loop QCD matrix elements is becoming well established, the virtual 
corrections form only part of a full NNLO calculation. They must be combined with the one-loop corrections 



to 7* — > 4 partons 1 14 - 17 , where one of the partons becomes collinear or soft, as well as tree-level processes 
7* — > 5 partons 1 1 1 - 1 3 1 with two soft or collinear partons in a way that allows all of the infrared singularities 
to cancel one another. This task is far from trivial, even though the factorization properties of both the one- 
loop, one-unresolved-parton contribution p7j-|72"| and the tree-level, two-unresolved-parton contributions J73] 
[76j have been studied. Early studies for the case of photon-plus-one-jet final states in electron-positron 
annihilation in J77|, [78|, which involves both double radiation and single radiation from one-loop graphs, 
indicate the feasibility of developing a numerical NNLO program implementing the experimental definition 
of jet observables and event-shape variables, and significant progress is anticipated in the near future. 

Note added: After this paper was first released, part of its results were confirmed in an independent 
calculation using the methods described in |4^,|^]]. In liep-ph/0207043, Moch, Uwer and Weinzierl obtain 



results for the full one-loop amplitude (3.14) and for the contributions to the two-loop amplitude ( |3.16 ) which 



are proportional to Np (i.e. the terms Dq and Eq), all in agreement with the results presented here. 
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A Weyl— van der Waerden spinor calculus 

The basic quantity is the two-spinor tpA or ip A and its complex conjugate tp A or ip A . Raising and lowering of 
indices is done with the antisymmetric tensor e, 



e AB =e AB =e AB = e AB = ^ ^ J j . (A.l) 
We define a antisymmetric spinorial "inner product" : 

(V'lV'a) = 1plAS BA 1p2B = 1plA^ A = -^t^2A = -0>2lpl), (A.2) 

and 

(1>iW = 1> lA 1>2- (A.3) 
Any momentum vector gets a bispinor representation by contraction with cr M : 

k AB = <* *„ = ( *° + *? k > + l t 2 ) , (A.4) 

AB AB v- y ki - ik 2 k - k 3 J ' K ' 

where er° is the unit matrix and cr. are the Pauli matrices. Since 



1 



we have 

For light-like vectors one can show that 
where 



k AB p AB = 2k-p. (A.6) 
k AB = k A k B, (A.7) 



u. — I (fci -i k 2 )/Vk ~k 3 . . . 



so that for light-like vectors we have 



2k-p^ (kp)(kp)* = \(kp)\ 2 . (A.9) 
The following relation is often useful: 

a» AB tf D = 25 A d 5 B D . (A.10) 
For massless spin-i particles the four-spinors can be expressed in two-spinors as follows: 



u+(p) = v_(p) = ( ^ ) ' 



-0) = v+(p) 







m+(<7) = v-(q) = ( 0, -ig^ ) , 

0_(g) - «+(«) = ( V 1 , ). (All) 
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The Dirac 7 matrices now become 



■ m 



r= \ ., AB * a , (A.12) 



so that, for example: 

u + (q)rv-(p) = q A a» AB p B . (A. 13) 

The general electroweak vertex for vector boson V coupling to two fermions is denoted by ieSijYy 1 * 2 , 
where i and j are the colour labels associated with the fermions fi and f2 respectively. The vertex contains 
left- and right-handed couplings, 

<v,hh = T v . „ ( 1 ~ 75 ^ ^ ,,r . i 1 + 

where for a photon, 
and for a Z-boson, 



~ j/ 1 - sin 2 6> we/l - -sine we/l 

£ ^= sin^cos^ Shh ' oose w Sf ^ (A ' 16) 

Here, ef represents the fractional electric charge, J3 the weak isospin and Oy/ the weak mixing angle. In the 
Weyl-van der Waerden notation, the vertex rjf'^ 1 ^ 2 becomes, 



" ~ iR\ f „<r AB ■ (A - 17) 



—iL>X, f„cr 

For the polarization vectors of outgoing gluons and photons we use the spinorial quantities 

6 iB (fc) = ^W' (A - 18) 

The gauge spinor b is arbitrary and can be chosen differently in each gauge-invariant expression. A suitable 
choice can often simplify the calculation. 
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